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Description 

BACKGROUND OF THE 




5 1. Field of the Invention. 

[0001 ] The present Invention relates generally to digital signal processing, and more particularly to signal processing 
for detecting the amplitude and frequency of a noisy, biased and sparsely-sampled sinusoidal signal. 

10 2. Background Art 

[0002] A known method of estimating the amplitude of a sampled sinusoidal signal is to compute a "root-mean square" 
(RMS) value, and then multiply the RMS value by the square root of two. This method, however, is inaccurate when 
the sinusoid is sparsely sampled and an estimate is obtained from a block of a small number of samples. 

15 

SUMMARY OF THE INVENTION 

[0003] In accordance with a basic aspect of the invention, there is provided a method of estimating the amplitude 
and frequency of a sinusoidal signal represented as sample values S 1 , S 2 , S 3> S 4 at respective spaced instants in time. 

20 The method includes the steps of: computing a triplet of differences x 1t x 2 , x 3 , where x 1 is a difference between S 2 and 
S 1t x 2 is a difference between S3 and S 2 , and x 3 is a difference between S 4 and S 3 ; and computing an estimate (B) of 
the amplitude, and an estimate (f) of the frequency, of the sinusoidal signal from the triplet of differences. Preferably 
an indication of the estimate of the amplitude orf requency is computed as a ratio of algebraic functions of the differences 
x 1? x 2 , x 3 . The amplitude, for example, is computed as the square root of a first ratio of algebraic functions, and the 

25 frequency is computed as an arc-cosine function of a second ratio of algebraic functions. Preferably the ratio is eval- 
uated by a division operation computed as a polynomial approximation, and the square root and the arc-cosine functions 
are also computed as polynomial approximations. 

[0004] In accordance with another aspect, the present invention provides a method of obtaining an indication of a 
characteristic of a signal represented as sample values (x 1f X2, x 3 , x n ) at respective spaced instants in time. The 

so method includes the steps computing a series of ratios of algebraic functions of triplets (x |} x M , Xj_ 2 ) of neighboring 
ones of the sample values, and computing an average of the ratios to produce the indication of the signal characteristic. 
Preferably any ratio having a denominator value that has a magnitude less than a threshold value is excluded from 
the average of the ratios, so that the average will not be corrupted by any inaccurate result that could be obtained 
when diving by a relatively small number. 

35 [0005] In accordance with a final aspect, the present invention provides a system for monitoring the health of a 
plurality of analog signals in an electronic circuit. The system includes an analog multiplexer for scanning the analog 
signals, an analog-to-digital converter for converting a selected analog signal to a series of digital values, a programmed 
data processor, and an output unit for indicating an error condition to an operator. The data processor is programmed 
to operate the analog multiplexer and to receive and process the series of digital values from the analog-to-digital 

40 converter. In particular, the data processor is programmed to determine an indication of amplitude and an indication 
of frequency of each of the analog signals by computing ratios of algebraic functions of at least four samples in the 
respective series of digital values for each of the analog signals, and to compare the indication of amplitude and the 
indication of frequency to limit values to indicate the error condition when one of the indications of amplitude and 
frequency falls outside of boundaries set by the limit values. 

45 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0006] Other objects and advantages of the present invention will become apparent upon reading the following de- 
tailed description and upon reference to the drawings, in which: 

so 

FIG. 1 is a block diagram showing a first method of using the invention, wherein the amplitude and frequency of 
a sinusoidal signal is estimated in order to provide a predetermined correction for parasitic capacitance in a circuit; 
FIG. 2 is a block diagram showing a second method of using the invention, wherein a multiplexer scans a number 
of signals in a system, and the amplitude and frequency of each selected signal is estimated and compared to 
55 predetermined limits in order to monitor the health of the system; 

FIG. 3 is a first portion of a flow diagram showing the present invention implemented for continuously estimating 
the amplitude and frequency of a sampled sinusoidal signal; 

FIG. 4 is a second portion of a flow diagram showing the present invention implemented for continuously estimating 
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the amplitude and freaug^y of a sampled sinusoidal signal; 

FIG. 5 is a third portic^^Bflow diagram showing the present invention implemWred for continuously estimating 
the amplitude and freqo^rrcy of a sampled sinusoidal signal; and 

FIG. 6 is block diagram showing the present invention implemented for obtaining an estimate of amplitude and 
5 frequency from a block of samples of a sinusoidal signal. 

[0007] While the invention will be described In connection with certain preferred embodiments, it Is not intended that 
the invention should be limited to these particular embodiments. On the contrary, the invention is intended to include 
all modifications, alternatives, and equivalent arrangements as may be included within the scope of this invention as 
10 defined by the appended claims. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

[0008] Turning now to the drawings, there is shown in FIG. 1 a block diagram illustrating a first method of using the 
*5 invention, wherein the amplitude and frequency of a sinusoidal signal is continuously estimated In order to provide a 
predetermined correction for parasitic capacitance in a circuit. In particular, FIG. 1 shows an angular rate sensor 10 
excited by a sinusoidal drive signal v d at a frequency f G , and producing an output signal v^ Indicating the level of 
excitation of the sensor, and a signal v r including angular rate information. The angular rate information, however, 
appears as a double-sideband suppressed-carrier signal about a suppressed carrier at the frequency f 0 . Such an 
20 angular rate sensor is described, for example, in Fersht et al. U.S. Patent 5,056,366 and Staudte U.S. Patent Re 
32,931 , herein incorporated by reference. 

[0009] As shown in FIG. 1 , the angular rate sensor 1 0 is excited by a driver circuit generally designated 1 1 so that 
the output signal v do has a predetermined amplitude. For this purpose, the signal v do is passed through a bandpass 
filter 1 2 to remove noise, and then the amplitude of the signal v^ is detected by an amplitude detector 1 3. An integrating 

25 comparator 14, having an associated input resistor 15 and a feedback capacitor 16, compares the amplitude of the 
signal v d0 to a predetermined reference V R , to produce an automatic gain control voltage V AGC . The automatic gain 
control voltage V AGC controls the gain of a variable-gain amplifier 17 to produce the exciting drive signal v d . 
[0010] To detect an angular rate signal a demodulator circuit 18 includes a balanced modulator 1 9 that is excited 
at the frequency f 0 . The balanced modulator 1 9 demodulates the signal v r to produce a base-band signal that is lowpass 

30 filtered by a lowpass filter 20 to produce the angular rate signal Q. 

[001 1] Although not pertinent to the present invention, further details regarding the angular rate sensor 10, the driver 
1 1 , and the demodulator 1 8 can be found in U.S. Patent 53 610 36 issued Nov. 1, 1994 , [Docket No. 92E083] by the 
present inventor Stanley A. White, and entitled "COMPLEX DIGITAL DEMODULATOR EMPLOYING CHEBYCHEWVP- 
PROXIMATION DERIVED SYNTHETIC SINUSOID GENERATION." 

35 [0012] In practice, the detected angular rate signal Q. is a function of the amplitude of the drive signal v d . Although 
the variation in the detected angular rate signal Q with respect to the amplitude of the drive signal v d is rather small, it 
is significant for some applications, and it is practical to remove the variation by employing the present invention to 
estimate the amplitude of the drive signal v d . For example, the variation in the detected angular rate signal Q with 
respect to the amplitude of the drive signal v d can be considered to be caused by a parasitic capacitance C p between 

40 nodes 1 9 and 20 in FIG. 1 . In this case, the error G in the detected angular rate signal Q, is approximately: 



E = KAf 0 (Equation 1) 

where K is a predetermined constant, A is the amplitude of the drive signal v d> and f c is the frequency of the drive 
signal. The constant K is on the order of 2TJC p Z in G i where C p is value of the parasitic capacitance, Zj n is the impedance 
at the node 22, and G is the conversion gain of the demodulator 18. 

[0013] In accordance with an aspect of the present invention, the error £= is determined by estimating at least the 
amplitude A of the drive signal v d , as further described below. The frequency f D could also be estimated, but since the 
frequency f 0 is relatively constant in comparison to the amplitude A of the drive signal v d , the frequency f 0 need not be 
estimated, and instead could be assumed to be a fixed, nominal value. The amplitude and/or frequency of the drive 
signal v d is estimated by sampfing the drive signal with an analog-to-digital converter 23, and processing the samples 
in an amplitude and frequency estimator 24. The amplitude and frequency estimator 24 will be further described below 
with reference to FIGS. 3 to 6. A correction signal generator 25 generates a correction signal e s based on the estimated 
amplitude and/or frequency, and an adder unit 26 subtracts the correction signal C s from the angular rate signal Q to 
produce a corrected angular rate signal Q\ 

[0014] Turning now to FIG. 2, there is shown a block diagram illustrating a second method of the invention, wherein 
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35 



an analog multiplexer 30j^s a number of signals S\ S 2 , S 3 , ... , S n from a 31, and the amplitude and 

frequency of each selectcflBal is estimated and compared to predetermined lim^^ order to monitor the health of 
the system. The system 3T^r example, is an inertial measurement unit including three angular rate sensors as de- 
scribed above with reference to FIG. 1 , each of the angular rate sensors have an associated driver circuit and detector 
s circuit, and each of the angular rate sensors sense the angular rate of rotation about a respective one of three orthogonal 
axes. I n this case the signal S 1 could be the drive signal v d for the rth one of the angular rate sensors, and the amplitude 
and/or frequency estimate for each of the drive signals could also be used for correction of the angular rate signal Q 
from each of the angular rate sensors. 

[0015] As shown in FIG. 2, the multiplexer 30 selects the signal S), and the selected signal is periodically sampled 

io by an analog-to-digital converter 32 at a sampling rate f s set by a sampling clock generator 33. The values of the 
samples are received at an input port 34 of a data processor 35. The data processor is interrupted at the sampling rate 
f s to receive the values of the samples from the analog-to-digital converter 32. An amplitude and frequency estimator 
program 36 controls the data processor 35 to receive a block of the sampled values, and to estimate the amplitude Al 
and the frequency fi of the sampled signal Si. A health monitoring program 37 controls the data processor 35 to compare 

15 the amplitude and frequency estimates to respective predetermine limits to monitor the health of the system. When 
the estimated amplitude or frequency is found to be outside the bounds of the predetermined limits, the data processor 
35 operates a display device 36 to display an error message to an operator (not shown). The health monitoring program 
37 then controls the data processor to change a selection signal from an output port 39 so that the analog multiplexer 
30 selects a different one of the signals S 1 , S 2 , S 3 S n . 

20 [0016] In a specific example, the signals S 1 , S 2 , S 3 , ... , S n each have a frequency of about 10 kilohertz, and the 
sampling rate is 42 kilohertz. The analog-to-digital converter is an Analog Devices part No. AD676 16-bit 100 kSPS 
sampling ADC which uses a switched-capacitor/charge redistribution architecture to achieve a 10 microsecond total 
conversion time. The data processor 35 is a Star Semiconductor SPROC programmable digital signal processing (DSP) 
integrated circuit which is designed for 24-brt arithmetic. In this system, the method of the present invention provides 

25 both amplitude and frequency estimates using a modest amount of computation and with errors on the order of parts 
per million. 

[001 7] Because the sampling rate f s is only about four (sometimes as low as three) times the frequency of the sinu- 
soidal signals S 1 , S 2 S 3 ... , S n , the sinusoidal signals are sparsely sampled. The amplitude and frequency estimation 
method of the present invention, however, provides an accurate estimate of the amplitude and frequency of such a 
30 sparsely sampled sinusoidal signal. The method of the present invention can also provide amplitude and frequency 
estimates based on as few as four samples of a sparsely sampled signal, in the presence of a DC bias. 
[001 8] Consider first the simplified case of a sinusoid without a DC bias. Assuming that the sinusoid has an amplitude 
A, a frequency o> 0 , and a phase angle and assuming that the sampling frequency f s has a period T=1/f s , then the n th 
signal sample x„ has a value of: 



40 



x n = A sin((0 o nT+(|)) (Equation 2) 

[0019] Since there are three unknowns in Equation 2, it should be possible to extract the amplitude A and a function 
of the frequency F=cosco 0 T from three samples of the sinusoidal signal. The two neighboring samples x^ and x n _.j are: 



X n+1 = A sin(co 0 nT+<l>)cos((o 0 T) + A cos(o 0 nT+<l>)sin(<D 0 T) (Equation 3) 

45 

X n-1 = A sin(a> 0 nT+<t>)cos(G> 0 T) - A cos(o 0 nT+<J>)sin(a> 0 T) (Equation 4) 

[0020] The function of frequency F is extracted from the three samples x n x„ +1 and x n . 1 by adding Equations 3 and 
4 together, and dividing the sum by Equation 2, to provide: 



x + x 

F= coso> 0 T= "* 1 °" 1 (Equation 5) 

[0021] The amplitude A is extracted by subtracting the product of Equation 3 and Equation 4 from the square of 
Equation 2, to provide: 
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X n - x n+1 x n-1 = A 2 Sin 2 (0 o T 




[0022] Applying the trigonometric identity 1 = cos 2 a + sin 2 a to Equation 5 gives: 



sin 2 o) 0 T 



10 [0023] Dividing Equation 6 by Equation 7 gives: 



(Equation 6) 



(Equation 7) 



15 



Y - Y Y 

^2 _ n A n+1 A n-1 



I 2 *n J 



1 



4 v(\i + \i) : 



(Equation 8) 



[0024] A direct implementation of Equation 8 to solve for the amplitude A therefore requires four multiplies, 3 adds, 
and 1 divide. 

20 [0025] Consider now the more complicated problem of finding the amplitude of a sinusoidal signal superimposed on 
a DC bias. The signal becomes: 



25 



30 



35 



40 



y n = DC + Bsin(o> 0 nT+<l>) 



(Equation 9) 



[0026] Since there are now four unknowns, at least four samples are needed to determine the amplitude B and 
function of frequency F = cosa> 0 T. The DC bias is eliminated by taking the differences x n between successive signal 
samples y n and y n _., as follows: 

x n = Vm' 3^n-i 

= B[sin(cD 0 nT+<|>) - sin([a> 0 (n-1)]7+<)>)] 
= B[sin(<o 0 nT+(|>) - sin (a> 0 /7T+<t>)cosa> 0 T + cos((0 o nT+<t»)sina) o 7] 
= B[sin(a> 0 nr+<t>) (1-cosa) 0 7) + cos(<o 0 n7+<|>)sina> 0 T| 
= e^/2-2cosa> 0 Tsin(a) 0 nr+e) 

where 



45 



8 = <|) + tan 



{ i-cos« 0 rj 



50 



(Equation 10) 

[0027] The DC biased sinusoid problem has been transformed to a mathematical form similar to the case without 
DC bias, except Equation 10 needs to be solved for B. But by comparison: 



55 



A= Bj2-2cos(i) 0 T 



(Equation 11) 



so that: 
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(Equation 12) 



s or: 



o2 A 2 *n' X n+i X n-1 



Solving Equations 6 and 7 for F gives: 



15 ir_ * n +t + X n-1 



Eliminating Ffrom Equation 14 gives: 



F= — p7 — - (Equation 14) 



_ _ (Equation 15) 



25 Factoring and re-grouping the terms in the denominator of Equation 1 5 gives an equivalent expression for the amplitude 
B: 

3 n o = o (Equation 16) 

[2 V(x n+1 + « 2x n) + < Vi + x n.i)l 

[0028] Turning now to FIGS. 3, 4, and 5, there is shown a fiow graph of a specific procedure implementing Equations 
14 and 1 6 above. This flow graph represents the solution procedure in a form similar to the representation for a digital 
35 filter, and indeed the elements in FIGS. 3, 4, and 5 could be individual computational hardware elements. Alternatively, 
the functions represented by the elements in FIGS. 3, 4, and 5 could be performed by a programmed digital data 
processor. 

[0029] As shown in FIG. 3, the sampled signal S, is delayed by one sample period T in a register 50 to produce the 
delayed sample S M . A subtracter 51 subtracts the delayed sample S M from the sample S, to produce a difference x,. 

40 The difference x, is delayed by a first sample period T in a second register 52 to provide a delayed difference x M , and 
is delayed for a second sample period T in a third register 53 to produce a doubly-delayed difference x,_ 2 . 
[0030] To compute the denominator of Equation 14, the delayed difference x M is scaled by a factor of two, as rep- 
resented by a box 54, which could left-shift the delayed difference by one binary place. To produce the numerator 
(num,) of Equation 14, an adder 55 adds the difference x, to the doubly-delayed difference x,_ 2 . 

45 [0031] To compute the denominator (DEN,) of Equation 16, a subtracter 56 computes the difference (den, - num,), 
and an adder 57 computes the sum (den, + num,). A multiplier 58 squares the difference from the subtractor 56, and 
a multiplier 59 multiplies the square by the sum from the adder 57 to produce the denominator DEN, of Equation 16. 
[0032] To compute the numerator (NUM,) of Equation 1 6, a multiplier 60 computes the square of the delayed differ- 
ence x M , and a multiplier 61 computes the cube of the delayed difference x M . A multiplier 62 computes the product 

so of the difference x, and the doubly-delayed difference Xj_ 2 . A subtractor 63 computes the difference between the square 
from the multiplier 60 and the product from the multiplier 62. A multiplier 64 computes the product of the difference 
from the subtractor 63 and the cube from the multiplier 61 . The product from the multiplier 64 is scaled by a factor of 
four as indicated in a box 65, which could represent a left-shift by two binary places, to produce the numerator NUM, 
of Equation 16. 

55 [0033] A divider 66 computes the quotient F, of Equation 1 3, and a divider 67 computes the quotient B, 2 of Equation 
15. The division performed by the divider units 66, 67 preferably is a polynomial approximation as derived in Appendix I . 
[0034] Continuing now in FIG. 4, the values of B, 2 are filtered in a digital filter 81 to determine average values B, 2 . 
The digital filter 81 , for example, is a non-recursive digital filter including a shift register storing the values B M 2 , B,_ 2 2 , 
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and B|_ 3 2 and computing tt^^im of Bj 2 and these three values and dividing the four, by a right-shift by two 

binary places, to provide a^vage value Bj 2 . The square root of this average valwiKomputed by a computational 
unit 82. The computational™ 82, for example, performs a square root operation by a polynomial approximation as 
derived in Appendix II. The result of the square root operation is an estimate B,of the amplitude of the sampled sinusoidal 
5 signal S. 

[0035] The division operation performed by the divider 76 in FIG. 3 may give an inaccurate result when the value of 
the denominator DENj is approximately zero. To prevent such inaccurate results from corrupting the amplitude estimate 
Bj, the inaccurate values of Bj 2 can be excluded so as not to be averaged by the digital filter 81 . For this purpose, a 
digital multiplexer 83 selects either the value Bj 2 from the division unit 67 in FIG. 3, or selects the previous average 

10 value computed by the digital filter 81 . The previous average value computed by the digital filter 81 is obtained from a 
register 84, that is clocked at the sampling frequency f s . To determine when the value B, 2 may be erroneous, the value 
DENj of the divisor for the division unit 67 is received by an absolute value unit 85 that determines the magnitude of 
the divisor value DEN,. A comparator 86 compares this magnitude to a certain fraction of the average value of the 
divisor value. A digital filter 87, which may be similar to the digital filter 81 , computes the average of the divisor value, 

15 and a shifter 88 right-shifts the average divisor value by N binary places to compute a predetermined fraction of the. 
average of the divisor value. The value of N, for example, is three. 

[0036] Turning now to FIG. 5, it is seen that the quotient F, of Equation 13 is processed in a similar fashion to the 
processing of the quotient Bj 2 as shown in FIG. 4, by a digital filter 91 , a functional unit 92, a multiplexer 93, a register 
94, an absolute value unit 95, a comparator 96, a digital filter 97, and a right-shifter 98, that are similar to the corre- 
20 sponding components 81 to 87 in FIG. 4. The functional unit 92 in FIG. 5, however, computes an arc-cosine function 
((1/2nT)acos(f ,)), The functional unit 92 preferably computes this arc-cosine function as a polynomial approximation, 
as derived in Appendix III. 

[0037] In practice, the estimates of amplitude Bj and frequency f, are relatively correlated between the neighboring 
ith and i+1th values. Therefore, it is possible to reduce the computational requirements of obtaining an estimate with 
25 only a slight decrease in performance by computing an amplitude estimate B and frequency estimate f from non- 
overlapping triplets of differences, rather than by computing the estimates from overlapping triplets of differences as 
illustrated in FIGs. 3, 4, and 5. In other words, although the computations illustrated in FIGs. 3, 4, and 5 can provide 
the best performance and provide a new estimate of the amplitude Bj and f ( at each sampling instant, the computational 
requirements for computing such estimates at each sampling instant are rather burdensome because ail of the corn- 
so putations shown in FIGS. 3, 4, and 5 must be performed over each sampling period T. For monitoring the hearth of 
signals in a system as illustrated in FIG. 2, for example, it is not necessary to compute a new estimate of amplitude or 
frequency for each sampling instant. Therefore, the computational requirements for monitoring the health of signals in 
a system as shown in FIG. 2 can be considerably relaxed by performing the computations illustrated in FIG. 3 over 
non-overlapping triplets of differences, and by performing the digital filtering and square root or arc-cosine function 
35 from the ratios computed from a number of non-overlapping triplets of differences. 

[0038] Turning now to FIG. 6, there is shown a flow diagram for the computations that are preferably performed for 
the system of FIG. 2 in the data processor 35. As illustrated in FIG. 6, an amplitude estimate B and a frequency estimate 
f are obtained from 1 6 non-overlapping triplets of differences x 1 , x^ x 3 ; x 4 , x 5 , x^ x 7 , x 8 , Xg; . . .; x^, x^, x^. The triplets 
of differences are the results of difference computations 101, 102, 103, ... 104, that are preferably performed by an 
40 interrupt program when the data processor 35 of FIG. 2 is interrupted at the sampling frequency f s . A quotient B, 2 and 
F, for each of the triplets of differences are computed by ratio computations 1 05, 1 06, 1 07 ... 1 08 that are performed 
at a rate of one-third of the sampling frequency f s . The 16 quotients Bj 2 and the 16 quotients f, are accumulated in 
digital filtering 109 that can be performed at various times over about 48 sampling periods T corresponding to the 
differences x, to x^. The average values of the quotient B 2 and F result from the digital filtering. A square root function 
45 110 computes the amplitude estimate B from the average value of the quotient B 2 , and the frequency estimate f is 
computed by an arc-cosine function 111 from the average value of the quotient F. 

[0039] Also shown in FIG. 6 is a variable Q associated with the amplitude estimate B, and a variable q associated 
with the frequency estimate f . The variables Q and q indicate, respectively, the number of the quotient values Bj 2 and 
Fj that were actually averaged together by the digital filter 109 instead of being excluded because their associated 

50 denominator value DEN, or denj, respectively, was relatively small. In contrast to FIGs. 4 and 5, where the exclusion 
condition was a certain fraction of an average value of DEN, or den,, respectively, the digital filter 1 09 preferably rejects 
quotients when the value of DEN or den falls below a respective predetermined threshold value. In the system of FIG. 
3, for example, the amplitude estimate B or the frequency estimate f should fall within certain predetermined limits, so 
that the threshold can be set at a certain fraction of the average value expected for any healthy signal. The threshold 

55 value is not critical and can be determined by simple trial and error for a particular system. 

[0040] Listed below is an example program that performs the functions shown in FIG. 6 and which corresponds to 
the amplitude and frequency estimator program 36 and the health monitoring program 37 of FIG. 2. In particular, in 
the listing below, there is shown a main program called "MONITOR" that operates the analog multiplexer 30 of FIG. 2 
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to scan among three signal^^Don outputting a new select signal to the multiplexe^Mfie program MONITOR waits 
30 microseconds for the at^^fco-digitat converter 32 to sample the new signal, T\mMe program MONITOR calls a 
routine called TESTSIG in ^^to test whether the amplitude and frequency of the selected signal are within prede- 
termined limits. If not, the program MONITOR displays an error code provided by the TESTSIG routine, and then the 
program MONITOR selects another signal for processing. 

[0041] The TESTSIG routine initializes a number of variables used by a 42 kHz interrupt routine corresponding to 
the amplitude and frequency estimator program 36 of FIG. 2. After initializing these variables, the TESTSIG routine 
enables a 42 kHz interrupt so that the amplitude and frequency estimator program processes the signal samples from 
the analog-to-digital converter 32 in FIG. 2. The routine TESTSIG polls a variable NUMC to determine when the am- 
plitude and frequency estimator program is finished. When NUMC is greater or equal to 16, then the amplitude and 
frequency estimator program 36 is finished. The TESTSIG routine then compares the variables Q and q to predeter- 
mined minimum limits QMIN and qmin, and also compares the amplitude estimate B and the frequency estimate f to 
respective predetermined minimum and maximum limits. An error indication is produced when any one of the variables 
Q, q, B, or f is outside of a boundary set by a respective one of the limits. 

[0042] The 42 kHz interrupt routine inputs one sample during each interrupt and classifies the sample as either 
caseO, easel , case2, or case3. CaseO corresponds to the very first sample S1 , and the other three cases correspond 
to one of the samples that permits the computation of either a first difference, a second difference, or a third difference, 
respectively, for each of the 16 triplets of differences. When a triplet of differences has been accumulated, then the 
computations depicted in FIG. 3 are performed for the triplet to compute a quotient B { 2 and Fj for the triplet. Also included 
in the program listing are examples of functions for computing the ratios, the square roots, and the arc-cosines by the 
polynomial approximations as derived in Appendices I, II, and III. 



MONITOR SELECT <- SELECT +1 

IF SELECT > 3 THEN SELECT <- 1 

OUTPUT SELECT 

WAIT 30 MICROSECONDS 

CALL TESTSIG 

IF ERROR > 0 THEN GOTO 10 

B (SELECT) <- SQRT(BS) 
*** NOTE: K = 1/2ITT 

f (SELECT) <- K*ACOS(F) 

RETURN 
10 DISPLAY(ERROR) 

RETURN 
TESTSIG ERROR <- 0 

CASE <- 0 

NUMC<-0 

BS<-0 

F<-0 

Q<-0 

q <- 0 

ENABLE 42 kHz INTERRUPT 
110 IF NUMC < 16 THEN GOTO 110 IF Q ^ QMIN THEN GOTO 120 ERROR <- 1 

RETURN 
120 BS<-BS/Q 

IF BS <> BSMAX THEN GOTO 130 ERROR <- 2 

RETURN 

130 IFBS>BSMIN THEN GOTO 140 ERROR <— 3 

RETURN 

140 IF q> qmin THEN GOTO 160 

ERRORS 4 

RETURN 
150 F<-F/q 

F < FMAX THEN GOTO 160 ERROR <- 5 

RETURN 

160 IF F > FMIN THEN RETURN ERROR <- 6 

RETURN 
42KH2INT INPUT SAMPLE 
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IF CASE > THEN GOTO 210 S1 <- SAMPLE 

CASE <- 1 
RETURN 

210 IF CASE > 1 • THEN GOTO 220 X1 <- S1 - SAMPLE 

5 S1 <- SAMPLE 

CASE<-2 
RETURN 

220 IF CASE > 2 THEN GOTO 230 X2 <- S1 - SAMPLE 

S1 *- SAMPLE 
10 CASE <- 3 

RETURN 
230 X3 <- S1 - SAMPLE 

S1 <r- SAMPLE 

CASE <- 1 
is num <- X1 + X3 

PRODX1X3<- X1 *X3 

NUMC<-NUMC + 1 

IF NUMC > 16 THEN GOTO 240 CLEAR INTERRUPT MASK 
GOTO 250 

DISABLE 42 KHZ INTERRUPT 
SQX2 <- X2 * X2 
CUBX2 <- SQX2*X2 
den 4- LEFT SHIFT X2 BY 1 
SUMnd<- num + den 
DEN «- den - num 
DEN <- DEN*DEN 
DEN <- DEN*SUMnd 
NUM <- SQX2 - PRODX1X3 
NUM <- NUMXUBX2 
30 LEFT SHIFT NUM BY 2 

IF den < denmin THEN GOTO 260 
F <- F + RATIO(num,den) 

q«-q + i 

260 IF DEN < DENMIN THEN GOTO 270 

35 BS <- BS + RATIO(NUM.DEN) 

Q<-Q + 1 
270 RETURN 



20 240 
250 



25 



40 



45 



50 



FUNCTION 
*** NOTE: 
*** B0 = 
*** B1 = 
*** B2 = 
*** B3 = 
*** B4 = 



55 



RATIO (T/U) 

0.8946079612 
-0.3153155446 
0.0547355413 
-0.0046808124 
0.0001578331 
RATIO <- U*B4 
RATIO <- RATIO + B3 
RATIO <- RATIO * U 
RATIO <r- RATIO + B2 
RATIO <- RATIO * U 
RATIO <- RATIO + B1 
RATIO <- RATIO * U 
RATIO <- RATIO + B0 
RATIO <- RATIO * T 
RETURN 



FUNCTION SQRT(V) 
*** NOTE: 
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15 



*** An — 
AU = 




*** A 1 — 




*** AO — 




*** A*3 — 

MO — 


U.UUDOO / O ID/ 


♦** A A — 

nf — 


"V.UVJUOUO/ 




oun i t — v D*+ 




own i < — oun i +■ do 




ouni t — ouni v 




CJORT i— QfiRT j_ R9 




SQRT<-SQRT* V 




SQRT<-SQRT+B1 




oun i t — own i v 






/ 


RPTI IRM 


pi IMPTIHM 




*** MHTP' 




*** rn — 


n nRnnnnniiQ 
u.uouuuuu 1 1 y 


0 1 = 


-u.DDoHDU^Hyy 


*** fsiy 

\jc. — 


nil nQ*;f;Annfc 


**« _ 


-u.uoouo l DO 




oLJVV < — VV W 




Ml»Uo <— oUW Oo 




MOUo < — MOUo + \jc. 




ACOS <— ACOS SQW 




ACOS <- ACOS + C1 




ACOS <- ACOS * W 




ACOS <- ACOS + CO 




RETURN 



30 

[0043] Finally, shown in Appendix IV is a computer program for evaluating the performance of various embodiments 
of the present invention. 

[0044] In view of the above, there has been described a method for estimating the amplitude and phase of a biased 
and sparsely-sampled sinusoidal signal. The method can give accurate results from a small number of samples, with 
35 only a modest amount of computation. Therefore, the method is easily implemented by a data processor programmed 
for estimating the amplitude and frequency of a number of analog signals scanned by an analog multiplexer and digitized 
by an analog-to-digital converter. 



40 



45 



50 
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APPENDIX I. 




MATHCAD (Trademark) PROGRAM 
CHEBYCHEV APPROXIMATION FOR DIVISION 
USING THE REMEZ- EXCHANGE ALGORITHM 
QRSDIVIDE . MCD 



We form the quotient x/y by evaluating 1/y and multiplying the 
result by x. The objective of this program is to calculate the 
°b M coefficients that provide the Chebychev solution to 1/y. 



INITIALIZE 

:= 6 Nominal value of the input 

variable . 

Tolerance on the input 
6 Y : = 1 .y aan variable - 



Y M := Y^ + 6Y = 8 Maximum anticipated 

input value, 

Y m» : = Y o«n " 6Y Y^ = 4 Minimum anticipated 

input value. 

B := 24 Number of bits of coefficient 

resolution, 

N := 5 Number of coefficients in 

data- fitting routine. (This 
is determined by trial and 
error. ) 

n : = 0 . .N Coefficient index. 

m := 0 . .N The y -value index. 

p := 0 ..N - 1 The weighting-coefficient 

(V) index. 

K : * 1000 Number of performance- 

evaluation points* 

k 0 ..K - 1 Performance -evaluation index. 



PRINTCOLWIDTH := 10 PRNPRECISION :« 10 

ZERO := 0 WRITEPRN(YA) :=*2ER0D 

m WRITE (COUNT) := 0 □ 

TYPE "PROCESS", COMMENT "WRITEPRN (YA) " & "WRITE (COUNT) M . THEN 

GOTO 125. 
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PERFORMANCE EVALUATION 



# 



AF 



The approximating funccion. 



ERROR := REF - AF 



The performance error. 



The plot of the error -extremal envelope bounded by € ' : 
c' := if [e = 0,max [i\ERROR\)l. Ml Another index i := 0 ..1 



This is the error -extremal envelope plot: 
EY k := if [k = K- 1,0, if[ERROR k . l ±ERROR k ,e',-€'\ ] 

This filter (MASK) is used to locate the frequencies 
corresponding to all extremal error values plus the end points: 

MASK k i= i£[k = 0,1, if [k * K - 1, 1. if [EY k , x = -EY k .l,0 ] ] ] 

This is the actual frequency- filtering operation. The only 
frequencies that get through are the end points plus those 
where error extrema are located. 

MATRIX ktC : « y k *MASK k MATRIX ktl : ' \ERROR k \ *MASK k 



Extract the frequencies where the N+l largest end-point values 

and extremal values are located, sort them, and write them into 
memory: 

MYE m .i reverse ( csor t {MATRIX, 1) ) a i SYE := (csort (MYE,0) ) 



|c| ■ 5.42 10' 5 WRITE PRN(YA) SYE <0> 
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max ((\ERH6R\) i^e\ = 1.25 10 



i-s 




A sanity check. 



floor 



-log [max [{\ERROR\)\ ] 



log (2) 



13 



Accuracy performance level 
measured in bits. 
Compute the peak error in 
ppm 



Now plot the functions and the 
error 



max [(\ERROR\)\ -10* = 66.7 



List the extremal frequencies 
(col 0) , and the extremal value 
magnitudes (col 1) : 



The rounded coefficients: 



The hex coefficients: 



SYE 



4 

14.3283283283 
S.24S2452452 
6.5025025025 
7 .5635635636 
8 



0.0000522733 
0.0000564838 
0.0000502018 
0.0000612455 
0.0000433821 
0.0000666976 



b = 



0.8946079612 
-0.3153155446 

0.054735S413 
-0.0046808124 

0.0001578331 



max 



SYE <X> \]\ _ 1 = 



min l\SYE< x> \]] 



0.5374449341 



2 B 'b 



eS0SO7h 
-S0b885h 
e0326h 
-132c3h 
aS8h 



END OF PROGRAM 



CHEBYCHEV APPROXIMATION FOR SQUARE ROOT USING THE REM EZ-EXCHANGE ALGORITHM QRSSQRT. 
MOD 

[0045] The objective of this program is to calculate the "a" coefficients that provide a Chebychev polynomial approx- 
imation to the square root of y. 

[0046] The Chebychev solution is the one that minimizes the maximum error of the approximation fit. There are no 
"surprise points" in the solution space. 
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+ 



INITIALIZE 

y • b 2 



Nominal value of the output 
variable. 



. v n , v Tolerance on the output 

or :- o.i^ variable. 



x ^ = 4.84 Maximum 
x oax ' 9 [ Y nom * 6 *j anticipated input value 



2 x ^ = 3.24 Minimum 

*iun [ r ooo * dy ] anticipated input value 



B :«24 Number of bits of 

coefficients resolution. 

N := 5 Number of coefficients in 

data-fitting routine. (This 
is determined by trial and 
error. ) 

n := 0 . .N Coefficient index. 

m := 0 ..N The y-value index. 

p := 0 ..N - 1 The weighting-coefficient 

("b") index. 

K := 1000 Number of performance - 

evaluation points. 

k : = 0 . . K - 1 Performance-evaluation index. 

PRINTCOLWIDTH 10 PRNPRECISION : = 10 

ZERO := 0 WRITEPRN(YA) : = ZERO □ 

m WRITE (COUNT) : = 0 O 

TYPE "PROCESS", COMMENT "WRITEPRN (XA) " & "WRITE (COUNT) " , THEN GOTO 
128. 
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SETUP 



ITERATION := READ (COUNT) WRITE (COUNT) 1 + ITERATION XA 

READPRN (XA) 

x H := if {iteration* o,a^ d ♦ [x^ - J^S.x^J 



The extreme values of the 
independent variable. 

The Remez matrix. 



- 1 



The performance - evaluation 
variable. 



R is the functional 
reference vector; REF is the 
performance -reference 
vector. 

v is the solution vector; € 
is the error measure 



A p := 2' B • fioor (0.5 + 2 B -v p ] 



The coefficient vector 
rounded to B bits. 
WRITEPRN ( ACOEFF ) a 



PERFORMANCE EVALUATION 



AF i = 



The approximating 
function. 



LP 



ERROR := REF - AF The performance error. 

The plot of the error-extremal envelope bounded by e': 
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€< := if [€ *Q,mmfrmmjT], |c|] 



J^p^ WOT ,_ M Another ind» i :-0 . .1 




This is the error- extremal envelope plot: 
EX k if [Jc - K - 1,0 if (ERKOR^* ERROR k ,s! , -c 7 ]] 



This filter (MASK) is used to locate the frequencies 
corresponding to all extremal error values, plus the end points: 

MASK* if [* « 0,1, if [Jc * JC - 1, 1 # if [EX^ * -£X*,1,0]J] 



This is the actual frequency- filtering operation. The only 
frequencies that get through are the end points plus those where 
extremal are located. 

MATRIX ktC : * x k • MASK k MATRIX ktl : = \ERROR£ • MASJC* 



Extract the frequencies where the N+l largest end-point values 

and extremal values are located, sort them, and write them into 
memory . 

MXE m ti reverse (csort (MATRIX, I) ) a i SXE (csortiMXE, 0) ) 

-6 <0> 
| € j , 1 #1 10 WRITEPRN(XA) ;= SXE 



—■ — — — - , . , ^ A sanity check, 

max [{\ERROR\)\ - |c| * 3.41 10 7 



Accuracy performance level 

flooi 



4 



log (2) 



19 measured in bits. 
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Now plot the fui^^ons and the Compute th^^ak error in 

error ppro: 

max[(|ER*0tf|)J -10 6 = 1.44 



List the extremal frequencies (col The rounded coefficients 

0) , and the extremal values (col are: 

1) : 



SXE 



3.24 
3.384144141 
3-7669269269 
4 .2682282282 
4.675035035 
4.84 



0. 0000010195 
0.0000011859 
0.0000009604 
0.0000013106 
0.0000007948 
0.0000014362 



a = 



0.5439969897 
0.5490743518 

-0.0689262748 
0.0068873167 

-0.0003057122 



""f fpg^ll l - ! , 0.80699 51557 
min ( [ iSXE^l ] ] 



The hex coefficients: 



2 B -a = 



8Jb4363ia 
Sc9023h 
-lla527i} 
lc35eA 
-1409J3 



END OF PROGRAM 



APPENDIX III. 



MATHCAD (Trademark) PROGRAM CHEBYCHEV APPROXIMATION TO AN ARCCOSINE FUNCTION 
QRSARCOSRMZ.MCD 

[0047] The objective of this program is to calculate the "a" coefficients that provide the Chebychev solution to the 
polynomial approximation of arccos(x). The application is a direct-reading frequency-deviation meter. We are going to 
estimate 5f/fo. 
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f, := 42000 
f 0 := 10000 
6f := 0.2-£ 0 

yL., : = cos 



INITIALIZE 




The sampling frequency. 

The nominal frequency. 

Maximum deviation of the 
variable to be estimated. 

Maximum value of the input 
variable Y^ = 0.3653410244. 



I min 



: - cos 



Minimum value of the input 
variable Y^ *= 0.222520934. 



B : = 24 



N := 4 



n 
m 
P 



0 
0 
0 



.N 
. N 

.N - 



1000 



k : = 0 . .K - 1 

PRINTCOLWIDTH := 10 
ZERO := 0 
m 

TYPE " PROCESS M , COMMENT 
THEN "GOTO 140". 



Number of bits of coefficient 
resolution . 

Number of coefficients in 
data- fitting routine. 

Coefficient index. 

The x- value index. 

The weighting-coefficient 
("a") index. 

Number of performance - 
evaluation points. 

Performance -evaluation index. 

PRNPRECISION := 10 
WRITEPRN(YA) := ZERCD 
WRITE (COUNT) : = 0D 
'WRITEPRN(YA) " & "WRITE (COUNT) rt . 



18 



. EP 0 638 811 B1 



SETUP 




ITERATION : = READ (COUNT) WRITE (COUNT) := 1 + ITERATION 

YA := READPRN (YA) 



Y m : = if {ITERATION * 0, * [Y m - J^J ^J, KA ffl ] 



The extremal values of the 
independent variable, y. 



The Remez matrix, 



variable. 



R is the functional reference 
vector; REF is the 
performance-reference vector. 



REF := - — 2_ -acos(y) - 1 



v := M" l -R € := v K 



a p :« 2~*-floor [0.5 * 2 B -v p J 



v is the solution vector; € is 
the error measure* 

The coefficient vector rounded 
to B bits. 

WRITEPRN (CCOEFF) := a 
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PERFORMANCE EVALUATION 




b p :» if[p * 0,0, a p ] 



ajt • = a* ♦ Pi> -v 2 ** 1 } The approximating 

«* \p u P y j function. 



ERROR :*= REF - AF The performance error. 

The plot of the error -extremal envelope bounded by e': 

c' := if [€ * Q,max[{\ERROR\)\ t \e\] Another index: i := 0 . 



EY k : * if [Jc * K - 1,0, if ERROR k ,*'. -e'] ] 



This filter (MASK) is used to locate the frequencies 
corresponding to all extremal error values: 

MASK k if [k = 0,1, if [k * K - 1,1, if [EY k . t * -EY k ,l,0] ] ] 



This is the actual filtering operation. The only frequencies 
that get through are those where error extrema are located. 

MATRIX^ := Y x • MASK^ MATRIX^ := |ERR0Rt| *MASK k 

Extract the frequencies where the N+l largest extremal values 
are located, sort them, write them into memory, and lisc the 
rounded coefficients: 

MYE ma := reverse (csort (MATRIX, 1))^ SYE : = 

(csort(MYE.O)} 
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|c[ = 5.02 10' 7 1 
max[(\ERROR\)\ - |«| = 1.58 10" 8 



floor 



-log [max {(\ERROR\) \ ] ] _ 2Q 
log(2) j " 



WRITEPRN (Y. 



9 



:= SYE 



<o> 



"Gut" insurance. (If you 
can't feel good about the 
solution, its probably 
wrong ! ) 

Accuracy performance level 
measured in bits. 



Now plot the functions and the 
error 



max [(\ERROR\)\ = 5.18 10* 7 



The extremal values: 



SYE 



-0.0830581871 
0.0817079274 
0.2294089799 
0.3294455494 
0.3653410244 



0.0000005136 
0.0000004898 
0.0000005165 
0.0000004871 
0,0000005183 



The rounded coefficients: 



3L - 



0.0500000119 
0.6684602499 
0.1108554006 
-0.0580515265 



max If \SYE <X> 1\ 



0.0639659227 



2 3 a = 



ccccdh 
-ab2036h 
-lc€105h 
-edc77h 



END OF PROGRAM 
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APPEMDIX IV. 



MATHCAD (Trademark) PROGRAM 
SIMPLE AND ROBUST QRS DRIVE -VOLTAGE 
AMPLITUDE AND FREQUENCY ESTIMATION 
QRSSINEAMPY.MCD 



SUMMARY 



This program introduces an algorithm that computes the peak 
amplitude of the sparsely sampled noisy sinusoid riding on a dc 
offset. The algorithm uses only adds, multiplies, and shifts 
operating on sets of 4 contiguous samples- The inherent 
computational error is on the order of 15 ppm. The algorithm 
also estimates frequency; the inherent computational error is 
within 2 or 3 hundredths of a Hertz. Naturally, high noise 
levels degrade the performance. 



SETUP AND DEFINITIONS 



K : = 


128 


k : = 


0 


. . K 


- 1 


Source -signal index 


L :« 


K - 1 


1 : = 


0 


. .L 


- 1 


Processed-signal index 


M : = 


L - 2 


m : = 


0 


. .M 


- 1 


Signal -processing index 






n : = 


0 


. .M 


- 1 


. . . and another one . 


*. 


42000 










Sampling frequency 


T : = 


1 


T = 


2. 


381 


io- s 


Sampling period 



I 

\ ' 5 2 A is the signal amplitude 

(the nominal value is Ao; 
we allow +/- 10% 
variation) 

A := A 0 -(1 + (-0.1 + rnd(0.2))) A = 1.937 

f 0 := 10000 Nominal signal frequency. 

f' 0 : = f 0 •(! + 0-2-{rnd(2) - 1)} Signal frequency (Hz). A 

+/- 20% variation is 
allowed about the 
nominal „ 
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Q := 2'ii'f' 0 Signal frequency 

(rad/sec) 

0 := rnd(2-n) * = 5.266 A completely arbitrary 

phase angle. 

:= A-sin(n-T-k + 0) QRS voltage sinusoidal - 

signal samples. 

?j k := 2- (rnd(2) - 1) Now we generate some 

noise ... 

DC := -1 + rnd(2) DC = -0.421 . and a DC offset 

(arbitrary) ... 

v : = u + tj + DC ... and add them to the 

signal in order to obtain 
the noisy and biased 
input voltage, v. 



SNR :« 10 -log 



k *K 



Signal-to-noise ratio {in 
dB) at the input: SNR = 
61. 



DEFINE SOME FUNCTIONAL APPROXIMATIONS 

We use Chebychev (minimax solutions) polynomial approximations to 
the square root, division, and arccosine operations. The 
polynomial coefficients have been rounded to 24 bits. These 
coefficients are imported from the SQRT.MCD, DIVIDE. MCD, and 



ARCOSRMZ.MCD programs. 



a := READPRN ( ACOEFF ) 



b := READPRN (BCOEFF) 



lasc(a) 



a ■ 



0.5439969897 
0.5490743518 
-0.0689262748 
0.0068873167 
-0.0003057122 



0.8946079612 
-0.3153155446 

0.0547355413 
-0.0046808124 

0.0001578331 



sgrt(x) : » Sa/x i 
sgn(y) : = if (y 2 0,1, -1) 
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div(x, y) : = s^A) x-2 




2- floor flSalMll f 2 -flooA l3liMl ] Y 
Z1 ° X \ log(2) \ Zb \\yY2 I log(2) J] 



e.g. , div(12,3) = 4.00V 

sqrt(4) = 2 sqrt(div(12,3)) = 2 



C : = READPRN (CCOEFF ) 



nfd = normalized frequency deviation. 



c = 



0.0500000119 
-0.6684602499 
-0.1108554006 
-0.0580515266 



Example \ 

iafd [cos [2n-£ 0 •(! + 0.1)- T] ] * 0.100001 



THE PROCEDURE 

Our first step is to eliminate the DC offset by taking the 
difference between successive samples of the input signal. The 
input signal to the differencing operator is v. the signal after 
the differencing operator is x: 

x i :=s v i+t - v, The difference signal 



The estimates of A and frequency are EoA : = J 2^ o Eof : = — —-a 

\ DEN 2ttT 

Define: 

DEN a :=[ - [x a + x^f ] • [2^ x - [x a ♦ x m . 2 } ) 

num* := x^ a + den^ := 2x tt+1 
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In order to avoTO division problems that arise if the 
denominators are too small, .we define the "weeding -out" 
functions, p and q: 

P a if[\DENj < 0.01- max [r\DEN\T), 0,1 ] 
g a if [\denj < 0,01-tnax \{\den\) J,0, 1 ] 

Now we compute the estimate of the square of A over several sets 
of samples r average those results, then take the square root: 

EoA :=sgr 



and compute the amplitude-estimation error in ppm: 
A 



We similarly estimate the frequency (the actual value is f'o) 
Eof : = £ 0 * |l * nfd J-JL ••aTvTHumTc^ J J 



.Eof - 9937.361 
f fi = 9937.135 



and determine the frequency-estimation error in Hertz: 

€ f f - Eof € f = -0.226 

O 



We have been operating with a signal-to-noise ratio of SNR = 
61.445 the phase angle is * = 5.266 the DC offset is DC = -0.421 
the frequency is off nominal by f - f = -62.865 and the 
amplitude is off by A - A = -0.063 o o. 

END OF PROGRAM 



Claims 

1 . A method of obtaining and using an indication of a characteristic of sinusoidal signal exciting an angular rate sensor, 
said method comprising the steps of: 
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a) exciting said an< 

b) digitizing said si 




rate sensor with said sinusoidal signal; 
fcal signal in an analog-to-digital converter, thereby 




icing a sequence of sample 



values S1 , S2, S3.^^^t respective spaced instants of time; 

c) computing in a computer operating pursuant to software a series of ratios of algebraic functions of triplets 
xi, xi-1 , xi-2 of neighboring ones of said sample values; 

d) computing in a computer operating pursuant to software an average of said ratios to produce said indication 
of said characteristic; 

e) applying said indication of said characteristic to a correction signal generator to produce a correction signal; 
and 

f) subtracting, in an adder unit, the correction signal from an angular rate signal generated by said angular 
rate sensor to produce a corrected angular rate signal. 

2. The method as claimed in claim 1 , wherein: 

the indication of the characteristic is an estimate (B) of the amplitude; 

the signal is a sinusoidal signal represented as sample values S1 , S2, S3, S4 at respective spaced instants 
of time; 

the number of sample values is four; and 

x1 is the difference between S2 and S1, x2 is the difference between S3 and S2, and x3 is the difference 
between S4 and S3. 

3. The method as claimed in claim 2, wherein said step of computing an estimate (B) of said amplitude is performed by 

i) computing a ratio of algebraic functions of said triplet of differences x1 , x2, and x3 ; and 

ii) computing said estimate (B) of said amplitude as a nonlinear function of said ratio of algebraic functions. 

4. The method as claimed in claim 3, wherein said ratio of algebraic functions includes terms equivalent to: 



5. The method as claimed in claim 4, wherein said nonlinear function includes a polynomial approximation of a square 
root function. 

6. The method as claimed in claim 2, wherein each of said steps a) and b) are performed by a data processor under 
program control. 

7. The method as claimed in claim 2, further comprising the step of computing a correction value from said estimate 
(B) of the amplitude of said sinusoidal signal for compensating a parasitic capacitance in an electrical circuit. 

8. The method as claimed in claim 1 , wherein: 

the indication of the characteristic is an estimate (F) of the frequency; 

the signal is a sinusoidal signal represented as sample values S1 , S2, S3, S4 at respective spaced instants 
of time; 

the number of sample values is four; and 

x1 is the difference between S2 and S1 , x2 is the difference between S3 and S2, and x3 is the difference 
between S4 and S3. 

9. The method as claimed in claim 8, wherein said step of computing an estimate (F) of said frequency is performed by 

i) computing a ratio of algebraic functions of said differences x1 , x2, and x3 ; and 

ii) computing said estimate (F) of said frequency as a nonlinear function of said ratio of algebraic functions. 

10. The method as claimed in claim 9, wherein said ratio of algebraic functions includes terms equivalent to: 



x n (*n ~ x n+1 x n-l) 
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11. The method as claimed in claim 10, wherein said nonlinear function includes a polynomial approximation of an 
arc-cosine function. 

12. The method as claimed in claim 8, wherein each of said steps a) and b) are performed by a data processor under 
program control. 

13. The method as claimed in claim 1 , wherein: 

a plurality of analog signals are produced at a like plurality of locations within an electronic circuit, one location 
for each respective signal; 

the plurality of signals are applied to an analog multiplexer which has a like plurality of inputs, one for each 

respective signal, and which further has an analog output and a selection control input; 

the analog output is connected to an input of an analog-to-digital converter which has a digital output; 

the digital output is connected to an input port of a programmed data processor which has a first output port 

connected to the selection control input of the analog multiplexer; 

the programmed data processor has a second output port which is connected to an output unit for indicating 

an error condition to an operator; and 

the programmed data processor is programmed to: 

1 ) output selection control signals to the analog multiplexer for conversion of each of the analog signals 
to a respective series of digital values by the analog-to-digital converter; 

2) input into the programmed data processor the respective series of digital values for each of the analog 
signals; 

3) determine an indication of amplitude and an indication of frequency of each of the analog signals by 
computing ratios of algebraic functions of at least four samples in the respective series of digital values 
for each of the analog signals; and 

4) compare the indication of amplitude and the indication of frequency to limit values to indicate the error 
condition through the second output port when one of the indications of amplitude and frequency falls 
outside a boundary set by one of the limit values. 

14. The method as claimed in claim 13, wherein said programmed data processor as programmed to compute said 
ratios of algebraic functions by computing a triplet of differences between successive values in each of said series 
of at least four samples, and by computing one of said ratios of algebraic functions indicative of amplitude from 
said triplet of differences, and computing one of said ratios of algebraic functions indicative of frequency from said 
triplet of differences. 

15. The method as claimed in claim 14, wherein said programmed data processor is programmed to average together 
a plurality of said ratios of algebraic functions to compute said indication of amplitude, and to average together a 
plurality of said ratios of algebraic functions to compute said indication of frequency. 

16. The method as claimed in claim 15, wherein said ratios of algebraic functions include algebraic functions which 
determine values for a denominator of each of said ratios of algebraic functions, and wherein said programmed 
data processor is programmed to exclude from said averages any ratios having a value for the denominator having 
a magnitude that is less than a limit value. 

17. The method as claimed in claim 1 , wherein each of said ratios has a denominator determined by a first one of said 
algebraic functions, and a numerator determined by a second one of said algebraic functions, and wherein said 
method further includes computing a value for said first one of said algebraic functions for each of said triplets, 
and excluding from said average any of said ratios for each of said triplets having a value for said first one of said 
algebraic functions that has a magnitude less than a threshold value. 

18. The method as claimed in claim 17, wherein said threshold value is constant for each of said triplets. 

19. The method as claimed in claim 17, wherein said threshold value is an average value of previously-computed 
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values for said first on^^iid algebraic functions 

20. The method as claim ec^raaim 1 7, wherein said, triplets are non-overlapping such at no two of said triplets both 
include any one of said sample values. 

21 . The method as claimed in claim 1 , wherein the three samples values xi, xi-1 , xi-2 in each of said triplets are 
computed from four sample values Sj, Sj-1, Sj-2, Sj-3 representing another signal, such that xi is a difference 
between Sj and Sj-1 , xi-1 is a difference between Sj-1 and Sj-2, and xi-2 is a difference between Sj-2 and Sj-3. 

22. The method as claimed in claim 1 , further comprising the steps of comparing said average value to a predetermined 
limit value, and when said predetermined limit value is exceeded, reporting an error condition to an operator. 

23. The method as claimed in claim 1 , wherein said ratio includes terms equivalent to: 

x n( x n~ x n+-\) x n-1 
4*„ 2 -(*„ + 1+*n-l) 2 

so that said ratio is an indication of amplitude of said sinusoideal signal. 

24. The method as claimed in claim 1 , wherein said ratio includes terms equivalent to: 

X n+1 + X n-1 



so that said ratio is an indication of frequency of said sinusoidal signal 

25. The method as claimed in claim 1, wherein the three samples values xi, xi-1, xi-2 in each of said triplets are 
computed from four sample values Sj, Sj-1, Sj-2, Sj-3, representing another signal, such that xi is a difference 
between Sj and Sj-1 , xi-1 is a difference between Sj-1 and Sj-2, and xi-2 is a difference between Sj-2 and Sj-3 and 
wherein said ratio includes terms equivalent to: 



[ 2 V(*n + 1+*n-l)l 2 [(2*nM*n + 1 +X n.l)] 



so that said ratio is an indication of amplitude of a sinusoidal signal, when said another signal represents said 
sinusoidal signal superimposed on a DC bias. 

26. The method as claimed in claim 1 , wherein each of steps a) and b) are performed by a data processor under 
program control. 



PatentansprUche 

1. Ein Verfahren zum Eriangen und Einsetzen einer Angabe bzw. eines Anzeichens fur eine Charakteristik eines 
sinusformigen Signals, das einen Winkelgeschwindigkeitssensor (angular rate sensor) anregt, wobei das Verfah- 
ren die folgenden Schritte aufweist: 

a) Erregen des Winkelgeschwindigkeitssensors mit dem sinusformigen (sinusoidal) Signal; 

b) Digitalisieren des sinusformigen Signals in einem analog-zu-digital Wandler, wodurch eine Sequenz von 
Abtast- bzw. Samplewerten S1, S2, S3....Sn an jeweiligen beabstandeten Zeitaugenblicken erzeugt wird; 

c) Berechnen in einem Computer, der gemaB einer Software betrieben wird, einer Serie von Verhaltnissen 
von algebraischen Funktionen von Dreiergruppen bzw. Triplets xi, xi-1, xi-2 von benachbarten Werten der 
Samplewerte; 

d) Berechnen in einem Computer, dergemaG Software betrieben wird, eines Durchschnitts der Verhaltnisse, 
urn eine Angabe der Charakteristik zu erzeugen; 
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e) Anlegen der An; 
zeugen; und | 

f) Subtrahieren, inc 




durch den Winkelgeschwindigkeitssensor, um ein korrigiertes Winkelgeschwindigkeitssignal zu erzeugen. 

2. Das Verfahren nach Anspruch 1 , wobei 

die Angabe der Charakteristik eine Schatzung (B) der Amplitude ist; 

das Signal ein sinusformiges Signal, dargestellt durch Samplewerte S1, S2, S3, S4 an jeweilig beabstandeten 
Augenblicken ist; 

die Anzahl der Samplewerte 4 ist; und x1 die Differenz zwischen S2 und S1 , x2 die Differenz zwischen S3 und 
S2, und x3 die Differenz zwischen S4 und S3 ist. 

3. Das Verfahren nach Anspruch 2, wobei der Schritt des Berechnens einer Schatzung (B) der Amplitude ausgefuhrt 
wird durch 

i) Berechnen eines Verhaltnisses der algebraischen Funktionen der Dreiergruppe der Differenzen x1 , x2 und 
x3; und durch 

ii) Berechnen der Schatzung (B) der Amplitude als eine nicht lineare Funktion des Verhaltnisses der algebrai- 
schen Funktionen. 

4. Das Verfahren nach Anspruch 3, wobei das Verhaltnis der algebraischen Funktionen Terme beinhaltet, die aqui- 
valent sind zu: 



5. Das Verfahren nach Anspruch 4, wobei die nichtlineare Funktion eine Polynomnaherung einer Quadratwurzel- 
funktion aufweist. 

6. Das Verfahren nach Anspruch 2, wobei jede der Schritte a) und b) durch einen programmgesteuerten Datenpro- 
zessor ausgefuhrt wird. 

7. Verfahren nach Anspruch 2, das weiterhin den Schritt des Berechnens eines Korrekturwertes der Schatzung (B) 
der Amplitude des sinusformigen Signals aufweist um eine parasitare bzw. Storkapazitat in einer elektrischen 
Schaltung zu kompensieren. 

8. Verfahren nach Anspruch 1 , wobei: 

die Angabe der Charakteristik eine Schatzung (5) der Frequenz ist; 

das Signal ein sinusformiges Signal ist, das als abgetastete Werte S1 , S2, S3, S4 an jeweiligen beabstandeten 
Augenblicken bzw. Zeitpunkten dargestellt ist; 
die Anzahl der Samplewerte 4 ist; und 

x1 die Differenz zwischen S2 und S1 , x2 die Differenz wischen S3 und S2, und x3 die Differenz zwischen S4 
und S3 ist. 

9. Verfahren nach Anspruch 8, wobei der Schritt des Berechnens einer Schatzung (F) der Frequenz ausgefuhrt wird 
durch 

i) Berechnen eines Verhaltnisses von algebraischen Funktionen der Differenzen x1 , x2 und x3; und 

ii) Berechnen der Schatzung (F) der Frequenz als eine nichtlineare Funktion des Verhaltnisses der algebrai- 
schen Funktionen. 

10. Verfahren nach Anspruch 9, wobei das Verhaltnis der algebraischen Funktionen Terme aufweist, die equivalent 
sind zu: 



x n(*n X n+1 x n-l) 



[ 2 V(^1^,l)f K2*>(*n + 1 + Vl)] 
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5 11. Verfahren nach Anspruch 10, wobei die nichtlineare Funktion eine Polynomnaherung einerArcuscosinusfunktion 
aufweist. 

12. Verfahren nach Anspruch 8, wobei jeder der Schritte a) und b) durch einen programmgesteuerten Datenprozessor 
ausgefuhrt wird. 

10 

13. Verfahren nach Anspruch 1 , wobei: 

eine Vielzahl von Analogsignalen an einer gleich groBen Vielzahl von Positionen innertialb eines elektroni- 
schen Schaltkreises erzeugt werden, jeweils eine Position fur jedes jeweilige Signal; 
is die Vielzahl der Signale an einem Analogmultiplexer angelegt werden, der eine gleich groBe Anzahl von Ein- 

gangen hat, jeweils einen fur jedes jeweilige Signal, und der weiterhin einen analogen Ausgang und einen 
Auswahlsteuereingang hat; 

der Analogausgang mit einem Eingang eines analog-zu-digital Wandlers verbunden ist, der einen digitalen 
Ausgang hat; 

20 die digitale Ausgabe verbunden ist mit einem Eingangsanschluss eines programmierten Datenprozessors, 

der einen ersten Ausgangsanschluss hat, der verbunden ist mit dem Auswahlsteuereingang des Analogmul- 
tiplexers; 

der programmierte Datenprozessor einen zweiten Ausgangsanschluss hat, der verbunden ist mit einer Aus- 
gabeeinheit zum Anzeigen eines Fehlerzustands an einen Betreiber; und 
25 der programmierte Datenprozessor programmiert ist zum: 

1 ) Ausgeben von Auswahlsteuersignalen zu dem Analogmultiplexer fur die Konvertierung eines jeden der 
Analogsignale in eine jeweilige Serle von digitalen Werten durch den analog-zu-digital Wandler; 

2) Eingeben der jeweiligen Serien von digitalen Werten in den programmierten Datenprozessor fur jedes 
30 der Analogsignale; 

3) Bestimmen einer Angabe der Amplitude und einer Eingabe der Frequenz fur jedes der Analogsignale 
durch Berechnen von Verhaltnissen von algebraischen Funktionen von zumindest vier Samples in den 
jeweiligen Serien von digitalen Werten fur jedes der Analogsignale; und 

4) Vergleichen der Angaben der Amplitude und der Angabe der Frequenz mit Grenzwerten, urn den Feh- 
35 lerzustand uber den zweiten Ausgabeanschluss anzuzeigen, wenn eine der Angaben der Amplitude und 

der Frequenz auBerhalb der Grenze, die durch die Grenzwerte gesetzt wird, liegt 

14. Verfahren nach Anspruch 13, wobei der programmierte Datenprozessor programmiert istzum Berechnen der Ver- 
haltnisse der algebraischen Funktionen durch Berechnung einer Dreiergruppe von Differenzen zwischen aufein- 

40 anderfolgenden Werten in jeder der Serie von zumindest vier Samplen und durch Berechnen eines der Verhaltnisse 

der algebraischen Funktionen anzeigend fur die Amplitude aus der Dreiergruppe der Differenzen, und Berechnen 
eines der Verhaltnisse der algebraischen Funktionen anzeigend fur die Frequenz aus der Dreiergruppe der Diffe- 
renzen. 

45 15. Verfahren nach Anspruch 14, wobei der programmierte Datenprozessor programmiert istzum Mitteln bzw. Durch- 
schnittsbilden einer Vielzahl der Verhaltnisse der algebraischen Funktionen, urn die Angabe der Amplitude zu 
berechnen, und zum Mitteln einer Vielzahl von den Verhaltnissen der algebraischen Funktionen, urn eine Angabe 
der Frequenz zu berechnen. 

50 16. Verfahren nach Anspruch 15, wobei die Verhaltnisse von algebraischen Funktionen algebraische Funktionen auf- 
weisen, die Werte fur einen Nenner fur jedes der Verhaltnisse von algebraischen Funktionen bestimmen, und 
wobei der programmierte Datenprozessor so programmiert ist, daB er von den Durchschnitten alle Verhaltnisse 
ausschlieBt, die einen Wert fur den Nenner bzw. Divisor haben, dessen Betrag geringer ist als ein Grenzwert. 

55 17. Verfahren nach Anspruch 1 , wobei jedes der Verhaltnisse einen Divisor hat, der durch eine erste der algebraischen 
Funktionen bestimmt ist, und einen Zahler hat, der durch eine zweite der algebraischen Funktionen bestimmt ist, 
und wobei das Verfahren weiterhin aufweist, dass ein Wert fur die erste der algebraischen Funktionen fur jede der 
Dreiergruppen berechnet wird und weiter von dem Durchschnitt all die Verhaltnisse fur jede der Dreiergruppen 
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i \^^fur die erste der algebraischen Funktionen habenM^e 



ausschlieBt, die einen ^^fur die erste der algebraischen Funktionen haberj^|^einen Betrag kieiner als ein 
Schwellenwert hat. 



18. Verfahren nach Anspruch 17, wobei der Schwellenwert fur jede Dreiergruppe konstant ist. 

5 

19. Verfahren nach Anspruch 1 7, wobei der Schwellenwert ein Durchschnittswert von zuvor berechneten Werten fur 
die erste der algebraischen Funktionen ist. 



20. Verfahren nach Anspruch 1 7, wobei die Dreiergruppen sich nicht uberlappen , so dass keine zwei der Dreiergruppen 
10 beide einen gleichen Abtastungswert beinhalten. 



21 . Verfahren nach Anspruch 1 , wobei die drei Samplewerte xi, xi-1 , xi-2 in jeder der Dreiergruppen aus vier Abtast- 
werten Sj, Sj-1 , Sj-2, Sj-3, die ein weiteres Signal darstellen, berechnet werden, so dass xi eine Differenz zwischen 
Sj und Sj-1 , xi-1 eine Differenz zwischen Sj-1 und Sj-2, und xi-2 eine Differenz zwischen Sj-2 und Sj-3 ist. 

15 

22. Verfahren nach Anspruch 1 , das weiterhin die Schritte des Vergleichens des Durchschnittswertes mit einem vor- 
bestimmten Grenzwert und Berichten eines Fehlerzustandes an einen Betreiber aufweist, wenn der vorbestimmte 
Grenzwert uberschritten wird. 



20 23. Verfahren nach Anspruch 1, wobei das Verhaltnis Terme aufweist, die aquivalent sind zu: 



4x*-(x n+1 +x /H ) 2 

25 

so dass das Verhaltnis eine Angabe uber die Amplitude des sinusformigen Signals ist. 
24. Verfahren nach Anspruch 1 , wobei das Verhaltnis Terme aufweist, die aquivalent sind zu: 

30 

X /741 + x n-1 



so dass das Verhaltnis eine Angabe uber die Frequenz des sinusformigen Signals ist. 

35 ^ 

25. Verfahren nach Anspruch 1 , wobei die drei Samplewerte xi, xi-1 , xi-2 in jeder der Dreiergruppen aus vier Sample- 
werten Sj, Sj-1 , Sj-2, Sj-3, die ein weiteres Signal darstellen, berechnet werden, so dass xi eine Differenz zwischen 
Sj und Sj-1, xi-1 eine Differenz zwischen Sj-1 und Sj-2, und xi-2 eine Differenz zwischen Sj-2 und Sj-3 ist, und 
wobei das Verhaltnis Terme aufweist, die aquivalent sind zu: 

40 

x n(*n~ x n+1 x n-l) 

[ 2x n-( x n + 1+ x n-l)l 2 [(2 x n H(X n+1 +X^ 1 )] 



45 so dass das Verhaltnis eine Angabe der Amplitude eines sinusformigen Signals ist, wenn das weitere Signal das 

sinusformige Signal uberiagert auf ein DC-Bias darstellt. 

26. Verfahren nach Anspruch 1 , wobei jeder der Schritte a) und b) durch einen programmgesteuerten Datenprozessor 
ausgefuhrt wird. 

50 

Revendications 

1 . Proc6d6 d'obtention et d'utilisation d'une indication d'une caract6ristique d'un signal sinusoidal excitant un capteur 
55 de vitesse angulaire, ledit proc6d6 comprenant les etepes consistant a : 

a) exciter ledit capteur de vitesse angulaire avec ledit signal sinusoidal ; 

b) num6riser ledit signal sinusoidal dans un convertisseur analogique/ num6rique, produisant ainsi une s£- 
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25 
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35 



quence de valeurs^^hantillonnage S1 , S2, S3, ... Sn k des instants temnfe espaces respectifs ; 

c) calculer dans u^^Bnateur fonctionnant conformement k un logiciel J^^pie de rapports de fonctions 
aigebriques de tripSRci, xi-1 , xi-2 de valeurs voisines desdites valeurs d'echantillonnage ; 

d) calculer dans un ordinateur fonctionnant conformement k un logiciel une moyenne desdits rapports pour 
produire ladite indication de ladite caracteristique ; 

e) appliquer ladite indication de ladite caracteristique k un g6n6rateur de signaux de correction pour produire 
un signal de correction ; et 

f) soustraire, dans une unite d'additionneur, le signal de correction d'un signal de vitesse angulaire g6n6r6 par 
ledit capteur de vitesse angulaire pour produire un signal de vitesse angulaire corrige. 

2. Procede selon la revendication 1 , dans lequel : 



indication de la caracteristique est une estimation (B) de I'amplitude ; 

le signal est un signal sinusoidal represents comme des valeurs d'echantillonnage S1, S2, S3, S4 k des 
*5 instants tempo re Is espaces respectifs ; 

le nombre de valeurs d'echantillonnage s'eieve a quatre ; et 
- x1 est la difference entre S2 et Si, x2 est la difference entre S3 et S2 et x3 est la difference entre S4 et S3. 

3. Procede selon la revendication 2, dans lequel ladite etape de calcul d'une estimation (B) de ladite amplitude est 
20 effectu6e en 

I) calculant un rapport de fonctions alg6briques dudit triplet de differences x1 , x2 et x3 ; et 

ii) calculant ladite estimation (B) de ladite amplitude comme une fonction non lineaire dudit rapport de fonctions 

aigebriques. 



4. Procede selon la revendication 3, dans lequel ledit rapport de fonctions alg6briques comprend des termes equi- 
valents k : 

x n( x n~ x n+1 x n-l) 



Wn+l+VlH 2 [(2X n WVl +X n-l)]' 



5. Proc6d6 selon la revendication 4, dans lequel ladite fonction non Iin6aire comprend une approximation polynomiale 
d'une fonction de racine carree. 

6. Precede selon la revendication 2, dans lequel chacune desdites etapes a) et b) est effectuee par un processeur 
de donnees sous le contr6le d'un programme. 

7. Procede selon la revendication 2, comprenant en outre retape de calcul d'une valeur de correction k partir de 
40 ladite estimation (B) de I'amplitude dudit signal sinusoidal pour compenser une capacite parasite dans un circuit 

eiectrique. 

8. Procede selon la revendication 1 , dans lequel : 

45 - indication de la caracteristique est une estimation (F) de la frequence ; 

le signal est un signal sinusoidal represente comme des valeurs d'echantillonnage S1 , S2, S3, S4 k des 

instants temporels espaces respectifs ; 

le nombre de valeurs d'echantillonnage s'eieve k quatre ; et 

x1 est la difference entre S2 et Si, x2 est la difference entre S3 et S2 et x3 est la difference entre S4 et S3. 

50 

9. Procede selon la revendication 8, dans lequel ladite etape de calcul d'une estimation (F) de ladite frequence est 
effectu6e en 

i) calculant un rapport de fonctions aigebriques desdites differences x1 , x2 et x3 ; et 
55 ii) calculant ladite estimation (F) de ladite frequence comme une fonction non lineaire dudit rapport de fonctions 

aigebriques. 

10. Procede selon la revendication 9, dans lequel ledit rapport de fonctions aigebriques comprend des termes 6qui- 
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valents&: ; cS^ifi 
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11. Procede selon la revendication 10, dans lequel ladite fonction non Iin6aire comprend une approximation polyno- 
miale d'une fonction d'arc cosinus. 

10 12. Procede selon la revendication 8, dans lequel chacune desdites etapes a) et b) est effectuee par un processeur 
de donnees sous le controle d'un programme. 

13. Procede selon la revendication 1, dans lequel : 

15 ' un © plurality de signaux anatogiques est produite k une pluralite identique d'emplacements k I'interieur d'un 

circuit electronique, avec un emplacement pour chaque signal respectif ; 
- la pluralite de signaux est appliqu6e k un multiplexeur analogique qui comporte une pluralite identique d'en- 
trees, une pour chaque signal respectif, et qui comporte en outre une sortie analogique et une entree de 
commande de selection ; 

20 - Ia sortie analogique est connectee k une entree d'un convertisseur analogique/ nurrterique qui comporte une 

sortie nurrterique ; 

la sortie nurrterique est connectee k un port d'entr6e d'un processeur de donnees programme dont un premier 
port de sortie est connects k I'entree de commande de selection du multiplexeur analogique ; 
le processeur de donnees programme comporte un deuxifcme port de sortie qui est connects k une unite de 
25 sortie pour indiquer une condition d'erreur k un op6rateur ; et 

le processeur de donnees programme est programme pour : 

1 ) deiivrer des signaux de commande de selection au multiplexeur analogique pour f aire convertir chacun 
des signaux analogiques en une serie respective de valeurs nurrteriques par le convertisseur analogique/ 

30 nurrterique; 

2) entrer dans le processeur de donnees programme la serie respective de valeurs nurrteriques pour 
chacun des signaux analogiques ; 

3) determiner une indication d'amplitude et une indication de frequence pour chacun des signaux analo- 
giques en calculant des rapports de fonctions algebriques d'au moins quatre echantillons dans la serie 

35 respective de valeurs nurrteriques pour chacun des signaux analogiques ; et 

4) comparer Indication d'amplitude et Vindication de frequence k des valeurs limites pour indiquer la 
condition d'erreur k travers le deuxfeme port de sortie quand I'une des indications d'amplitude et de fre- 
quence tombe en dehors d'une borne definie par I'une des valeurs limites. 

*o 14. Precede selon la revendication 13, dans lequel ledit processeur de donnees programme est programme pour 
calculer lesdits rapports de fonctions alg6briques en calculant un triplet de differences entre des valeurs succes- 
sives de chacun de la s6rie d'au moins quatre echantillons, et en calculant run desdits rapports de fonctions 
algebriques indiquant ['amplitude k partir dud it triplet de differences et en calculant I'un desdits rapports de fonctions 
algebriques indiquant la frequence k partir dudit triplet de differences. 



15. Proc6d6 selon la revendication 14, dans lequel ledit processeur de donn6es programme est programme pour 
etablir la moyenne d'une pluralite desdits rapports de fonctions algebriques pour calculer ladite indication d'am- 
plitude, et pour etablir la moyenne d'une pluralite desdits rapports de fonctions algebriques pour calculer ladite 
indication de frequence. 



16. Proc6de selon la revendication 15, dans lequel lesdits rapports de fonctions algebriques comprennent des fonc- 
tions algebriques qui determinent des valeurs pour un denominateur de chacun desdits rapports de fonctions 
algebriques, et dans lequel ledit processeur de donnees programme est programme pour exclure desdites moyen- 
nes tous les rapports comportant une valeur pour le denominateur ayant une grandeur qui est interieure k une 

55 valeur limite. 

17. Procede selon la revendication 1 , dans lequel chacun desdits rapports comporte un denominateur determine par 
une premidre desdites fonctions algebriques, et un numerateur determine par une deuxieme desdites fonctions 
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algebriques, et dans leougl ledit procede comprend en outre le calcul d'une v^^fcour ladite premiere desdites 
fonctions algebriques poumacun desdits triplets, et {'exclusion de ladite moyer^^tout rapport desdits rapports 
pour chacun desdits triplets ayant une valeur pour ladite premiere fonction desdites fonctions algebriques qui a 
une grandeur inferieure a une valeur de seuil. 

18. Precede selon la revendication 17, dans lequel ladite valeur de seuil est constante pour chacun desdits triplets. 

19. Procede selon ja revendication 1 7, dans lequel ladite valeur de seuil est une valeur moyenne de valeurs calculees 
auparavant pour ladite premiere fonction desdites fonctions algebriques. 

20. Procede selon la revendication 17, dans lequel lesdits triplets ne se recouvrent pas de sorte que deux desdits 
triplets n'incluent jamais a la fois Tune quelconque desdites valeurs d'echantillonnage. 

21 . Procede selon la revendication 1 , dans lequel les trois valeurs d'echantillonnage xi, xi-1 , xi-2 dans chacun desdits 
triplets sont calculees a partir de quatre valeurs d'echantillonnage Sj, Sj-1 , Sj-2, Sj-3 qui represented un autre 
signal, de sorte que xi est une difference entre Sj et Sj-1 , xi-1 est une difference entre Sj-1 et Sj-2 et xi-2 est une 
difference entre Sj-2 et Sj-3. 

22. Procede selon la revendication 1 , comprenant en outre les etapes consistant a comparer ladite valeur moyenne 
a une valeur limite predeterminee, et quand ladite valeur limite predeterminee est depassee, signaler une condition 
d'erreura un operateur. 

23. ProcedS selon la revendication 1 , dans lequel ledit rapport comprend des termes equivalents a : 

2, 2 . 
4x^(x n+1 +x fv1 ) 2 

de sorte que ledit rapport est une indication d'amplitude dudit signal sinusoidal. 

24. Procede selon la revendication 1 , dans lequel ledit rapport comprend des termes equivalents a : 

x m-1 + X n-1 
x n 

de sorte que ledit rapport est une indication de frequence dudit signal sinusoidal. 

25. Procede selon la revendication 1 , dans lequel les trois valeurs d'echantillonnage xi, xi-1 , xi-2 dans chacun desdits 
triplets sont calculees a partir de quatre valeurs d'echantillonnage Sj, Sj-1 , Sj-2, Sj-3 qui represented un autre 
signal, de sorte que xi est une difference entre Sj et Sj-1 , xi-1 est une difference entre Sj-1 et Sj-2 et xi-2 est une 
difference entre Sj-2 et Sj-3 ; et dans lequel le rapport comprend des termes equivalents a : 

x n( x n~ x n+1 x rvl) 

PV^n+l+VlM 2 KSXnWx^+X^)] 

de sorte que ledit rapport est une indication d'amplitude d'un signal sinusoidal quand ledit autre signal represente 
ledit signal sinusoidal superpose a une polarisation en courant continu. 

26. Procede selon la revendication 1 , dans lequel chacune des etapes a) et b) est effectuee par un processeur de 
donnees sous le controle d'un programme. 
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